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EXPERIMENTALLY DETERMINED LOCAL FLOW PROPERTIES AND 
DRAG COEFFICIENTS FOR A FAMILY OF BLUNT BODIES 
AT MACH NUMBERS FROM 2.49 TO 4.63 

By Robert L. Stallings, Jr. 

Langley Research Center 

SUMMARY 

Drag coefficients and local flow properties were experimentally determined for a 
family of blunt bodies at Mach numbers from 2.49 to 4.63 at a Reynolds number, based on 
afterbody diameter, of 1.88 x 10®. The family consisted of bodies of revolution having 
variable nose and shoulder radii (r n and r c , respectively) and cylindrical afterbodies 
7.5 inches (191 mm) in diameter (d). The geometry of the 18 models tested ranged from 
a hemisphere-cylinder to a flat-face cylinder. 

The Mach number effect on nondimensional pressure and velocity distributions of 
the hemispherical model decreased rapidly with increasing Mach number. These distri- 
butions at Mach number 4.63 were essentially the same as previously published results 
for Mach numbers up to 11.4. Increasing the nose bluntness also decreased the Mach 
number effect on the pressure and velocity distributions. There was no effect of Mach 
number on these distributions for the zero-shoulder -radius models for -^p> 0.707. Drag 
coefficients determined from integrated pressures over the nose of the hemispherical 
model and by assuming the base pressure coefficient to correspond to -i/m^ (where 
Moo is the free-stream Mach number) were in good agreement with previously published 
data for a sphere. These results were also in good agreement with drag coefficients 
determined from modified Newtonian theory. A reduction in drag coefficients as indicated 
by both experiment and theory occurred for a decrease in bluntness obtained by an 
increase in shoulder radius. The maximum velocity gradient for all models occurred 
either at or slightly downstream of the point of tangency of the nose and shoulder arcs. 
Stagnation-point velocity gradients determined from measured pressures were in good 
agreement throughout the range of variables of this investigation with theoretical esti- 
mates based on Traugott's method and measured shock-standoff distances. 

A comprehensive presentation of these data in figure form is included for suffi- 
ciently small intervals of nose and shoulder radii to enable the pressure distributions, 
velocity distributions, stagnation-point velocity gradients, shock -standoff distances, or 
drag coefficients to be determined — either directly or by interpolation — for any body 
of the general shape described. Since all these variables indicated only very small 
Mach number effects at the higher test Mach numbers, these results should be applica- 
ble at a much higher range of Mach number than that of this investigation. 



INTRODUCTION 


The advantages of blunt nose shapes for reducing the convective aerodynamic 
heating at the forward stagnation point of bodies moving at hypersonic flight speeds have 
been well established within the past decade. Unfortunately, the governing partial- 
differential equations for the subsonic flow field between the bow shock wave and nose of 
such bodies are of the elliptic type and, as yet, no exact analytical solutions are available. 
This condition has led to what is generally referred to as the "blunt body problem." 

Some gains have been made in recent years toward the solution of this problem by various 
numerical and approximate methods - for example, see references 1 to 6; however, as 
discussed in reference 7, large discrepancies can exist between the pressure distributions 
determined by the different methods. Until the problems associated with these discrep- 
ancies are resolved, experimental investigations are required to determine the local flow 
properties on all but the more basic nose shapes. 

A blunt nose shape that has received considerable attention in the past and that has 
been used on numerous reentry configurations (e.g., Mercury, Gemini, and Apollo) con- 
sists of a hemispherical segment nose with a shoulder region having a circular cross 
section. Experimentally obtained stagnation-point velocity gradients for bodies of this 
type with a shoulder radius of zero have been reported in reference 8 and with a limited 
range of shoulder radius, in reference 9. Although the stagnation-point velocity gradi- 
ents are extremely important for determining the stagnation -point heating, it is well 
known that for extremely blunt bodies the maximum heating occurs off the stagnation 
point. Therefore, in order to assess accurately the level of heating over such a family 
of bodies, detailed pressure and velocity distributions must be determined over the com- 
plete nose. Pressure distributions for such bodies have been reported in references 8, 

10, and 11; however, the models used in these investigations provided only a limited range 
of geometrical variables and, being small, were limited in the amount of instrumentation. 

The present investigation was therefore initiated to determine the effect of nose and 
shoulder radii on the local pressure and velocity distributions for a family of 18 blunt 
bodies. The models, which were 7.5 inches (191 mm) in diameter, ranged from a hemi- 
spherical nose to a flat -face cylinder at intervals of nose and shoulder radii sufficiently 
small to enable the results from this investigation to be applied - either directly or by 
interpolation - to any shape of the general type. The tests were conducted through a 
range of Mach number from 2.49 to 4.63. The flow properties presented and discussed 
for the complete range of geometrical variables consist of pressure distributions, veloc- 
ity distributions, stagnation -point velocity gradients, and shock- standoff distances. Drag 
coefficients obtained by integrating the local pressures are also discussed. Limited 
comparisons are made with approximate theories. 
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SYMBOLS 


a* free-stream sonic velocity 

A area 

Cj) drag coefficient 

Cp pressure coefficient 

Cp jin ax maximum pressure coefficient, based on total pressure behind normal shock 
D drag force 

d afterbody diameter 

M Mach number 

p pressure 

r^ afterbody radius 

r c shoulder (corner) radius 

r n nose radius 

s surface length, measured from forward stagnation point (see fig. 1) 

nose surface length from forward stagnation point to point of tangency of nose 

and shoulder arcs (see fig. 1) 

S 2 nose surface length from forward stagnation point to shoulder -afterbody 

juncture (see fig. 1) 

T temperature 

u velocity 

u f stagnation -point velocity gradient, du/ds 

0 angle between normal to model surface and axis of symmetry (see fig. 1) 

0\ value of 0 at sj 

y ratio of specific heats, 1.4 for air 
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shock- standoff distance 


p density 

Subscripts: 

co free-stream conditions 

l local static conditions at outer edge of boundary layer 

t free-stream stagnation conditions 

t,2 stagnation conditions behind normal shock at free-stream Mach number 

b pertaining to model base 

h pertaining to hemisphere 

APPARATUS AND TEST CONDITIONS 

The investigation was conducted in the high Mach number test section of the Langley 
Unitary Plan wind tunnel described in reference 12. This variable -pressure, continuous 
flow tunnel has an asymmetric sliding-block nozzle that permits a continuous variation 
in the test-section Mach number from 2.30 to 4.63. The deviations in Mach number in 
the entire 4- by 4-foot (1.22- by 1.22-meter) test section for the test Mach numbers are 


as follows: 

For Moo = 2.49 ±0.018 

For Moo =3.51 ±0.048 

For Moo = 4.06 ±0.061 

For Moo = 4.63 ±0.045 


The effects of nonuniform Mach number were minimized in this investigation by 
testing all models at essentially the same location in the test section. The effect of flow 
angularity associated with this nonuniform Mach number was also minimized by an adjust- 
ment of the models for each test point relative to the free-stream velocity vector. This 
step was accomplished by monitoring pressure differentials from the stagnation point to 
orifice locations equidistant and diametrically opposite the stagnation point and adjusting 
the model in both angle of attack and angle of yaw until these pressure differences 
equalized. 
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The pressure measurements were obtained for the model at an angle of attack of 0°. 
The free-stream stagnation temperatures at the test Mach numbers were as follows: 


Moo 

°R 

°K 

2.49 

610 

339 

3.51 

610 

339 

4.06 

635 

353 

4.63 

635 

353 


The test Reynolds number, based on afterbody diameter, was 1.88 x 10®. 

MODELS, INSTRUMENTATION, AND ACCURACY 

The general shape of the axisymmetrical models (see fig. 1) consisted of a 
hemispherical segment nose (of radius r n ) faired into a circular-arc shoulder (of 
radius r c ) which faired into a cylindrical afterbody (of diameter d). A total of 
18 models were tested and they had geometries ranging from a flat-face cylinder = 0, 
to a hemisphere =0.5, ^1 = 0.5^. Values of r c /d and r n /d for each 
model are shown in the table presented in figure 1. Also included in this table are values 
of sj/d and S 2 /d, where sj is the value of s at the point of tangency of the nose 
and shoulder arcs and S2 is the value of s at the shoulder -afterbody juncture. The 
afterbody for all models consisted of a cylindrical section 4 inches (102 mm) long and 
7.5 inches (191 mm) in diameter. The model instrumentation consisted of approximately 
80 pressure orifices of 0.050-inch (1.27-mm) inside diameter. Locations of these pres- 
sure orifices for a typical model are shown on the sketch in figure 1. Photographs of each 
model tested are shown in figure 2, and a typical model installation in the test section is 
shown in figure 3. 

In the forward stagnation region of blunt bodies, the pressure magnitudes are gen- 
erally quite large; however, the pressure gradients can be very small. Such a combina- 
tion makes it extremely difficult, if not impossible, to measure the magnitude of the pres- 
sure decrease with surface length within this region with an absolute-pressure gage to 
the precision required for accurately determining local velocity gradients. This problem 
was minimized in this investigation by using a sensitive differential-pressure gage, full- 
scale deflection of 1 psi (6895 N/m2), to measure the pressure differential from the for- 
ward stagnation point to a select number of locations. The magnitudes of the pressures at 
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Figure 3 - Typical model installation in test section. L-65-7350 

these locations were obtained, by subtracting the pressure differentials from the 
stagnation -point pressure which was measured with a precision mercury manometer. 

The pressures at all remaining locations were measured with an absolute transducer 
having a full-scale deflection of 10 psi (68950 N/m.2). Both the 1- and 10-psi (6895- 
and 68950-N/m.2) transducers were used in conjunction with a multichannel scanning 
system so that only a total of four transducers were required. The output from each 
electrical transducer was recorded with a digital self-balancing potentiometer. 

The tunnel free-stream static and total pressures were measured with precision 
mercury manometers. The accuracy of the precision mercury manometers is within 
0.5 psf (23.94 N/m.2); therefore, the accuracy of the pressure measuring system is lim- 
ited to that of the electrical transducers. The accuracy of the electrical transducers is 
within 1 percent of full-scale deflection, which corresponds to a pressure increment of 
1.44 and 14.4 psf (69 and 690 N/m.2) for the 1- and 10-psi (6895- and 68 950-N/m.2) 
gages, respectively. 

RESULTS AND DISCUSSION 
Pressure Distributions 

The effect of Mach number on the hemisphere pressure distributions is shown in 
figure 4 for the test range of Mach number from 2.49 to 4.63. The local measured pres- 
sures have been normalized by the measured pressure at s/d = 0. The effect of Mach 
number, as expected, consists of a decrease in the magnitude of the normalized pressure 
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Figure 4.- Comparison of pressure distributions obtained on hemispherical model with previously 

published results. 

distributions; the extent of this effect decreased with increasing Mach number. Also 
shown in this figure are experimental data from references 13 and 14 for Mach numbers 
up to 11.4. The data obtained from reference 13, indicated by the hatched region, cover 
a range of Mach number from 1.8 to 9.6. The Mach number 1.8 data are at the upper 
bound of the hatched region and the M = 9.6 data are near or at the lower bound of this 
region. There is only a very slight difference between the M = 4.63 data of the present 
tests and the higher Mach number data of both references 13 and 14. 

Presented in figure 5 is the effect of nose radius on the normalized pressure dis- 
tributions for the zero-corner-radius models. The two extreme geometries are the flat- 
face cylinder (tH and the hemisphere = 0.500^. The dashed parts of the curves 

for 1.000, 1.933, and «> are extrapolations from the pressure measurements on 

these models at the last instrumentation station to the pressure corresponding to sonic 
velocity, which for these models would occur at s = S 2 - For the other three models, the 
sonic point location occurred at s < S 2 * 
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Figure 5.- Effect of nose radius on pressure distributions of 
zero-corner- radius models. M® = 4.63. 


The effect of shoulder radius on the pressure distributions for constant values of 
r n jd is shown in figure 6. The test -point symbol notation corresponds to the same val- 
ues of r c /d for all parts of figure 6; however, the values of S 2 /d are a function of 
r n /d and therefore may have different values for the same symbol. The short-dashed 
curves for = 0 in figures 6(a), 6(b), and 6(c) are the same extrapolations shown in 
figure 5. The long-dashed curves labeled s^Jd represent the location on the pressure- 
distribution curves of the point of tangency of the nose and shoulder arcs for each model. 
The data represented by the circular symbols = 0.5 'j were obtained on the hemisphe- 
rical model and are therefore independent of r n jd. The curve for ^p = 0.5 also repre- 
sents one of the limiting cases for each value of r n /d. The other limiting case is the 
pressure distributions for the zero-corner -radius models, which are a function of r n /d. 
With decreasing values of r n /d, these distributions approach the distribution for the 
hemispherical model which reduces the extent of the corner-radius effect. 


A modified Newtonian pressure distribution was evaluated from the following 
expression for = 0: 


= cos2[S iL) + 

;.2 dr n/ 


sin2(§ <L 
l d r n 


10 









and is presented for each value of r n /d in figure 6. This theoretical distribution should 
also apply to those models having > 0 where < ^p since within this range of s/d 
the local surface slopes are the same. The limitation of the Newtonian theory for pre- 
dicting either the pressure magnitudes or the pressure distributions for 
when ^p < 0.5 is clearly indicated in figures 6(a), 6(b), and 6(c). As expected, for those 
models having = 0.707 (fig. 6(d)), much better agreement is shown with the Newtonian 
theory because the local surface slope at the shoulder juncture for these models is 
approximately the same as that at the sonic point for the hemispherical model. It is 
interesting to note that for each value of r n /d the pressure-distribution curves inter- 
sect each other at a common location (other than at the stagnation point) for values of 
0-2 = -^p = 0.5. The value of s/d at which this intersection occurs appears to increase 
slightly with decreasing values of r n /d. 

The effect of decreasing Mach number from 4.63 to 2.49 on the normalized pressure 
distributions is shown in figure 7 for selected values of r c /d at each value of r n /d. A 



Figure 7.- Effect of Mach number on pressure distributions. 
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comparison of these distributions, together with the Mach number effect shown for the 
hemispherical model in figure 4, clearly indicates that the Mach number effect rapidly 
diminishes with increasing nose bluntness (with decreasing shoulder radius). Further- 
more, since the results for the hemispherical model at M^ = 4.63 were in good agree- 
ment with higher Mach number results (from refs. 13 and 14), the data included in this 
paper for all models should be applicable at much higher Mach numbers so long as real- 
gas effects can be neglected. 


Drag Coefficients 

Drag coefficients were determined for each model through the test range of Mach 
number from integrated pressure measurements. Since the base of the models was not 
instrumented, the drag increment resulting from base pressure was estimated from the 
following often-used empirical expression: 

C D,h - Cp,b - - (2) 

With this relationship being assumed for Cp) } t>) a comparison of Cd for the hemisphe- 
rical model with existing experimental flight data for a sphere is shown in figure 8. 

Good agreement is shown between the present data and previously published data (refs. 15 
to 18). Also shown are the drag coefficients calculated for a base pressure equal to free- 
stream static pressure, as implied by Newtonian theory. These drag coefficients (solid 
circular test points) are considerably lower than the flight data at the lower test Mach 



Figure 8.- Variation with Mach number of drag coefficients for hemispherical model. 
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numbers although the magnitude of the disagreement decreases with increasing Mach 
number. The drag coefficients evaluated from the modified Newtonian theory 


Cpjinax 
= 2 


- c p,b 


( 3 ) 


are based on the standard Newtonian assumption 
C D b = — i-ir. Use of the -1/M 2 relation for 

HIT 4 ' 


Cp^b = 0 as well as on the expression 
Cp b, as shown in figure 8, gives much 


better agreement between experiment and theory at the lower Mach numbers. At the 

higher Mach numbers, the values of C n k = — i-s- approach zero and therefore the drag 

V ’ M^ 


coefficients approach the standard Newtonian values. 


The effects of shoulder and nose radii on the drag coefficients for all test Mach 
numbers are presented in figure 9. The effect of shoulder radius for a given nose radius 
can be seen by observing the variation of a given test-point symbol with r c /d whereas 
the effect of nose radius is shown by the different symbols at a constant shoulder radius. 



Figure 9.- Variation of drag coefficients with model nose geometry. 
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Theoretical values, determined from modified Newtonian theory, were evaluated from the 
following equations as derived in appendix A: 


£n <00 
r b 


C D “ 2^] C Pj max ^ ” cos ^ e l) + 2Cp,max f " ^ sin 0 l( cos ^ 0 l + 2 ] 


. 1 r rc 

+ 2 u P,max[— 


COS^0J + 


M. 


(4) 


and for ~ = », 
r u 


- ^p, max * 1 


^ +#C. 


3 '-'p, max 


£c/i _ £c 


4c 

+ 2 u P,max 




(5) 


For Moo = 2.49 (fig. 9(a)), the drag coefficients decrease with either decreasing nose 
radius or increasing shoulder radius because of a reduction in the extent of bluntness, 
however, the effect of nose radius decreases with an increase in shoulder radius. 

Although this trend is shown by both experiment and theory, the decrease in the effect 
of nose radius on Cp> with an increase in shoulder radius occurs much more rapidly 
for the experimental data. In fact, for ^- = 0.2, the experimentally obtained drag coef- 

Cl 2* 

ficients remained essentially constant for 0.707 ^ = °° whereas the theoretical drag 

coefficients decrease by approximately 9 percent as a result of this range of reduction 
in T n jd. 

For the flat -face cylinder model = oj, the experimental drag coefficient 

is less than the theoretical Newtonian value by approximately 8 percent. This disagree- 
ment results from the fact that the physical flow expands to sonic velocity at the shoulder 
with the corresponding reduction in local pressure (see fig. 6(a)) whereas the Newtonian 
theory predicts a constant pressure across the face. For = 0, the agreement between 

theory and experiment improves with decreasing r n /d so that for ^ 0.707 the dis- 
agreement is within the data accuracy. Intuitively, the agreement between theory and 
experiment would be expected to improve with increasing values of r c /d at a constant 
value of r n /d since the body geometry approaches the hemispherical shape for which 
the Newtonian theory is known to give good results. As shown in figure 9(a), this 
improvement occurs for values of > 0.2; however, for values of < 0.2, the theo- 
retical and experimental curves diverge with increasing r c /d for 0.707 «>. 

2* 

The maximum disagreement is approximately 16 percent and occurred for = 00 and 

^•=0.1 or 0.2. C 

d 



The same general trends in the variation of drag coefficient with model geometry 
as shown for M = 2.49 are shown in figures 9(b), 9(c), and 9(d) for Mach numbers 3.51, 
4.06, and 4.63, respectively. The results presented in figure 9(d) are replotted in fig- 
ure 10 in order to show more clearly the effect of nose radius on the drag coefficients. 
Note that the abscissa scale in figure 10 is the inverse of the parameter r n( /d used in 
the key of figure 9. This change was made in order to eliminate the broken scale that 
would be required to show the results obtained at = *>. As previously discussed, the 

variation of Cj ) with r n jd is very small for ^ = 0.2; however, these data, as pre- 
sented in figure 10, actually indicate a slight increase in Cj) with increasing d/r n 
for — - ^ 1.0. This trend, which is contrary to what might be expected (since increasing 
djr n decreases the extent of bluntness), is also indicated for values of ^ = 0.1 and 0.3. 


Velocity Distributions 

Velocity distributions for the hemispherical model, nondimensionalized by sonic 
velocity, are presented in figure 11 for the test Mach number range. Also shown in this 
figure are distributions obtained by previous investigators (refs. 13 and 14 and 19 to 22) 
for Mach numbers from 1.97 to 11.4. 

The referenced data for Mach num- 
bers from 1.97 to 6.8 as shown by the 
faired curves are in good agreement 
with the present results for both the 

Theory Exp. r c /d 
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Figure 10.- Effect of nose radius on drag 
coefficients. = 4.63. 


Figure 11.- Effect of Mach number on velocity distributions for 
hemispherical model. 
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velocity distributions and the magnitude of the Mach number effect on these distributions. 
The velocity distributions shown for Moo =9.6 (ref. 13) and Moo = 11.4 (ref. 14) are 
slightly less than those shown for either M^ = 6.8 (ref. 19) or Moo = 4.63 of the pres- 
ent tests. This apparent discrepancy of approximately 5 percent is probably within the 
accuracy of the referenced results for the higher Mach number. 

The nondimensional velocity distributions obtained on the zero-corner -radius 
models are shown in figure 12 for M^ = 4.63. The test -point symbols in this figure 
correspond to the same geometrical variables used for the pressure results presented 
in figure 5. In order to improve the accuracy of the velocity distributions for small val- 
ues of s/d, these as well as all subsequent velocity distributions were determined from 
averaged pressures around the face of the models. This averaging process could result 
in slight differences for velocities computed from the pressures presented in figure 5 
which are only for one instrumented ray. For all values checked, the deviation of the 
averaged velocity from those along any particular ray was well within the instrumentation 
accuracy. 



Figure 12.- Effect of nose radius on velocity distributions of zero-corner- 
radius models. Moo = 4.63. 
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The effect of shoulder radius on the nondimensional velocity distributions for con- 
stant values of r n jd is presented in figure 13. The dashed curves labeled sj/d again 
represent the location of the point of tangency of the nose and shoulder arcs for each 
value of r c /d as shown previously in the figures of pressure data. Although no appar- 
ent unique characteristic of the pressure-distribution curves occurred at sj/d, an 
inflection occurred in the velocity-distribution curves either at or slightly downstream 
of this location for ^ > 0 and all values of r n /d. This inflection corresponds to the 
maximum velocity gradient, the magnitude of which decreases with increasing r c /d at 
constant values of r n /d. As the values of r n /d are increased, the velocity-distribution 
curves for the two limiting models (with shoulder radii = 0 and 0.5^ tend to move 
toward each other and therefore reduce the overall effect of r c /d. This tendency is 
similar to that shown by the results for the pressure distributions. 

The effect of Mach number on the nondimensional velocity distributions is shown 

in figure 14 for selected geometrical shapes, as well as in figure 11 for the limiting 

hemispherical model. The extent of the Mach number effect decreases with increasing 

bluntness similar to that shown by the pressure -distribution results. The velocity- 

distribution results for the zero-corner -radius models with 0.707 show no Mach 

d 



Figure 13.- Effect of shoulder radius on velocity distributions for constant nose radius. M„ - 4.63. 
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number effect whereas those for the models with ^ = 0.707 show an increase with 
increasing Mach number. This increase is due to the fact that the sonic point for 

0.707 is fixed at the shoulder, whereas the sonic point for ^ = 0.707 has moved 
forward of this juncture so that the Mach number effect is essentially the same as for a 
hemisphere. 


Stagnation -Point Velocity Gradients 

Stagnation-point velocity gradients were determined graphically from large-scale 
plots of the averaged velocity measurements in the stagnation region of the models. Only 
those velocities evaluated from the pressures measured with the 1-psi (6895-N/m2) 
transducers were used. Typical velocity distributions and the linear fairings for the 
stagnation region are shown in figure 15. 

A correlation parameter derived from modified Newtonian theory and the incom- 
pressible Bernoulli equation was used to correlate the stagnation -point velocity gradients 
for the hemispherical model with those obtained from references 19, 20, and 22 to 24. 
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s, in. s, in. s, in. 

(a) Model 11; r c /d = 0 ; r n /d = 1.000. 



Figure 15.- Typical velocity distributions in stagnation region of models. 


From the modified Newtonian theory, the pressure distribution can be expressed as 


= cos^0 + -^°°- sin^0 (6) 

p t,2 p t,2 


and from the incompressible Bernoulli equation, which is valid within the model stagna- 
tion region, the velocities can be written in terms of pressures as follows: 



( 7 ) 
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Combining equations (6) and (7) and simplifying yields 


u = sin g 


2 ^t> 2 (i 


p t,2 \ P t, 2 


1/2 


Since for small values of 9, 


sin g « g 


( 8 ) 


the following expression is obtained: 


u * 9 


' 2p t ,2 L 


p t,2 \ p t,2 


1/2 


Differentiating equation (9) with respect to s, where 


, _ du dg _ 1 du 
dg ds “ r n dg 


the stagnation-point velocity gradient for air can be expressed as 

* 58.6 


r^u 

n 




1/2 


( 9 ) 


( 10 ) 


r n u ' 


The parameter on the left-hand side of equation (10) was evaluated for the stagnation- 
point velocity gradients obtained on the hemispherical model and is presented in figure 16 
as a function of Mach number. 

Also shown in figure 16 are 
results obtained from refer- 
ences 19, 20, and 22 to 24 for 
Mach numbers up to 6.8. The 
results from the present inves- 
tigation fall within the data 
scatter of the previously 
obtained results and all 
results are within ±7 percent 
of the constant value predicted 

by equation (10). Figure 16.- Variation with Mach number of stagnation-point velocity-gradient 

parameter for hemispherical model. 
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Since Newtonian theory has limited application to the other models of this investi- 
gation, the use of equation (10) was restricted to the hemisphere. In an attempt to mini- 
mize the Mach number effect on stagnation -point velocity gradients for the other models, 
the velocity gradients were nondimensionalized by the values obtained from the hemi- 
spherical model at each test Mach number. This parameter was used inasmuch as the 
stagnation-point velocity gradients for all models decreased monotonically with increasing 
Mach number. These nondimensional velocity gradients are presented in figure 17 for 
all models investigated through the test range of Mach number. Also shown in figure 17 

are data obtained from reference 8 for values of = 0 and ^ = °° and 1.0. These 

d d 

results are approximately 6 percent greater than those of the present tests. The values 
of u'jj used for nondimensionalizing the reference 8 data were determined from modi- 
fied Newtonian theory. However, as shown in figure 16, the measured values of u' for 
the hemispherical model for the present tests are approximately 7 percent greater than 
those obtained by the Newtonian theory and this discrepancy is believed to be the major 
source of the difference between the two sets of u/u'jj data. The theoretical values 



Figure 17.- Variation of stagnation-point velocity gradients with model geometry and Mach 
number. (Flagged symbols determined from eg. (11); solid symbols from ref. 8.) 
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of u/u'k shown in figure 17 were determined from the following equation as derived in 
appendix B: 


u' _ 6 h 1 + (Wr n ) 
u 'h 6 l + (8/r n ) h 


( 11 ) 


where the values of 6 were determined from schlieren photographs and are presented in 
nondimensional form in figure 18. The values of 6/d shown in figure 18 for ^ = 0.5 
are from previously obtained results for a 9 -inch -diameter spherical model tested in 
the Langley Unitary Plan wind tunnel since schlieren photographs were not obtained for 
the hemispherical model of the present investigation. Also shown in figure 18 are data 
from references 25 and 26. 



Figure 18.- Variation ot shock-standoff distance with model geometry and Mach number. 

(Solid symbols from ref. 25; flagged symbols from ref. 26.) 

In general, the values of u'/u'h determined from equation (11) are in good agree- 
ment with the experimental values. (See fig. 17.) It should be emphasized that u'h, 
as used in figure 17 and in all subsequent figures containing the parameter 
corresponds to a hemisphere having a radius equal to the afterbody radius of the general 
shape depicted in figure 1 and not equal to the nose radius of this general shape. 
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The nondimensional velocity gradients shown in figure 17 at a Mach number of 4.06 
are replotted in figure 19 to show the effect of shoulder radius more clearly. The param- 
eter u'/u'h is plotted in figure 19 as a function of r c /d for various values of r n /d. 
The advantage of this form of presentation is that for a given Mach number the velocity 
gradients for this complete family of blunt bodies, ranging from the zero-corner -radius 
models = 0^ to the hemispherical model = 0.5^, can be shown in one figure for all 
values of r n /d. Furthermore, the results in figure 19 should apply to Mach numbers 
higher than 4.06 because the nondimensional velocity distributions (see fig. 11 as well as 
fig. 14) are essentially invariant with further increases in Mach number. 

The nondimensional velocity gradients shown in figure 17 at a Mach number of 4.06 
are replotted in figure 20 to show the effect of nose radius more clearly. The param- 
eter u'/u'h is plotted as a function of nose radius for various values of r c /d. In order 
to retain the data for the flat -face model without using a broken scale, the velocity gra- 
dients are plotted as a function of d/r n rather than the normally used parameter r n /d. 
Also shown in figure 20 are the theoretical estimates from modified Newtonian theory 
(dashed line) and from equation (11) (flagged symbols). The expression for the modified 
Newtonian theory for this family of bodies reduces to 


u' _ ( rn )h .Id. 
u’ h r n “ 2 r n 


( 12 ) 


as can be seen from equation (10). The limited application of the Newtonian theory for 
predicting stagnation-point velocity gradient for bodies of this type is clearly indicated 
in figure 20; however, the estimates based on equation (11) are in good agreement with 
the data throughout the range of geometrical variables. 



Figure 19.- Effect of shoulder radius on 
stagnation-point velocity gradients. 

Mco = 4.06. 



Figure 20.- Effect of nose radius on stagnation- point 
velocity gradients. IWx> = 4.06. (Flagged symbols 
denote theoretical values from eq. (ID.) 
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SUMMARY OF RESULTS 

Drag coefficients and local flow properties were experimentally determined for a 
family of blunt bodies at Mach numbers from 2.49 to 4.63. The family consisted of bodies 
of revolution having variable nose and shoulder radii (r n and r c , respectively) and cylin- 
drical afterbodies 7.5 inches (191 mm) in diameter (d). The model geometry ranged 
from a hemisphere-cylinder ^^p =0.5, ^p = 0.5) to a flat-face cylinder ^p = », ^p = 0^. 
The results are summarized as follows: 

1. The Mach number effect on both the nondimensional pressure distributions and 
the velocity distributions of the hemispherical model decreased rapidly with increasing 
Mach number through the Mach number range of this investigation. These distributions 
at a Mach number of 4.63 were essentially the same as previously published results for 
Mach numbers up to 11.4. 

2. The Mach number effect on the pressure and velocity distributions decreased 

with increasing the extent of nose bluntness. There was no effect of Mach number on 

these distributions for the zero-corner-radius models for > 0.707. 

d 

3. Drag coefficients determined from integrated pressures over the nose of the 
hemispherical model and by assuming the base pressure coefficient to correspond to 
-l/Moo^ (where Moo is the free-stream Mach number) were in good agreement with 
previously published results for a sphere. These results were also in good agreement 
with drag coefficients determined from pressure distributions evaluated over the model 
nose by modified Newtonian theory and the use of the aforementioned expression for the 
base pressure coefficient. 

4. A reduction in drag coefficients as indicated by both experiment and theory 
occurred for a decrease in bluntness obtained by an increase in shoulder radius. The 
maximum disagreement between experiment and theory was approximately 16 percent 
and occurred for ^p = °° and ^p = 0. 1 or 0.2. There was essentially no variation in 
the drag coefficients with Mach number at the higher test Mach numbers. 

5. The maximum velocity gradient for all models occurred either at or slightly 
downstream of the point of tangency of the nose and shoulder arcs. 

6. Stagnation-point velocity gradients determined from measured pressures were 
in very good agreement throughout the range of variables of this investigation with theo- 
retical estimates based on Traugott's method and measured shock-standoff distances. 

7. A comprehensive presentation of these data in figure form is included for suffi- 
ciently small intervals of nose and shoulder radii to enable the pressure distributions, 
velocity distributions, stagnation-point velocity gradients, or shock-standoff distances to 
be determined — either directly or by interpolation — for any body of the general shape 
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described. Since all these variables indicated only very small Mach number effects at 
the higher test Mach numbers, these results should be applicable to a much higher range 
of Mach number than that of this investigation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 11, 1967, 

124-07-02-44-23. 
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APPENDIX A 


NEWTONIAN DRAG COEFFICIENTS 

For the family of axisymmetric bodies defined by the geometrical parameters shown 
in the following sketch: 


*c 




the differential drag resulting from a pressure p^ acting on a differential surface area 
dA can be expressed as 

dD = p^ cos 9 dA (Al) 

The total drag would then become 




Pi 


cos 9 dA 


-i 


Nose 


p. cos 9 dA + 


I, 


Base 


P l cos 9 dA 


(A2) 


Since the value of 9 is constant (180°) for the base of this family of bodies, the base 
drag component can be easily integrated if p^ is assumed to be constant. With this 
assumption, equation (A2) becomes 


Dm i 


Nose 


cos 9 dA - p^A^ 


(A3) 


where Aj, is the base area. 

In order to present equation (A3) in the form of pressure coefficients, the 
expression 



cos 9 dA = 0 


is used in conjunction with equation (A3) to obtain 

D = ^Nose " P ~ )COS 9dA -(^- P ~) A b 


(A4) 


(A5) 
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APPENDIX A - Continued 


In the form of drag coefficient, equation (A5) becomes 

Cd = Tr~ \ Cn cos 9 dA - C D b (A6) 

A b J Nose P P ’ 

From modified Newtonian theory, the pressure coefficient can be expressed as 

Cp = Cp,max cos 2 0 (A7) 

which when substituted into equation (A6) gives 

C D = ~ f cos 3 0 dA - C p b (A8) 

A b '-'Nose 

The differential area dA for the nose region can be expressed in terms of 9 as 
follows: 

For 0 S 0 g flj, 


dA = 27rr n 2 sin 0 d0 


(A9a) 


and for 6i = 9 = 


dA = 27 rr c (r b - r c + r c sin 9)d9 


(A9b) 


Combining equations (A8) and (A9) leads to the following integral expression for the drag 
coefficient: 



Integrating equation (A10) and rearranging terms yields the following closed-form 
equation: 



cos^l) + 2Cp >max 



£c\ 2 

r b/L 3 





+ «• Cp^maxfjr cos ^^l _ ^p,b 



(All) 
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APPENDIX A - Concluded 


P 

This equation is of indeterminate form for — = however, by taking the appropriate 

r b 

limiting process, it can easily be reduced to the following form: 


CD - Cp >max (l 




^1 


r b) 3 ' r£ 


■c^ + l r 

+ 2 c P,max(] 


8 


^ - C 


P,b 



(A12) 
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APPENDIX B 


NONDIMEN SIONAL STAGNATION- POINT VELOCITY GRADIENTS 


Using the first Belotserkovskii approximation, Traugott (ref. 5) obtained simple 
analytical relations for the stagnation region of blunt bodies between shock and body cur- 
vature, shock -standoff distance, and velocity gradient. These relations enable an approx- 
imation of stagnation-point velocity gradients to be made from optical information, which 
is generally much easier to obtain than pressure data. The equation from reference 5 
relating stagnation-point velocity gradient with shock -standoff distance is 


where 


1 

6 

Uoo 

1 + (6/r n ) 


vr(Moo,y) 


(Bl) 


ir(Moo,y) = 


2 + (y - 1)M 0 
(y + 1)M m 2 


i - r 


2(M 00 2 - l) 
[2 + (y - l)Mo 


p yMpp 2 - 2(y - 1) 
(y + 1) 2 M co 2 


il/(y-l) 


and 


TtM^y) = 


2(y + 1) + 2y(y + l)Moo 2 - (y + l) 2 M 0O 2 ((y + l) 2 M 0O 2 /[4yM oo 2 - 2(y - l)j] 


4y(M 0O 2 - l) 


Since both n and T are independent of model geometry, they can be eliminated 
from equation (Bl) by nondimensionalizing the stagnation-point velocity gradient by the 
gradient for some reference shape. A logical shape to use is the hemisphere since both 
the stagnation-point velocity gradient and shock-standoff distance are well defined through 
a wide range of Mach number. Solving equation (Bl) for u' and nondimensionalizing by 
u’ h gives 

u r _ 5h 1 + (S/r n ) 

u' h ~ a i + (a/r n ) h ^ 

which is applicable for any value of M w or y. 
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